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Abstract 



Fusion hierarchies of A-D-E face models are constructed. The 
fused critical D, E and elliptic D models yield new solutions of the Yang- 
Baxter equations with bond variables on the edges of faces in addition 
to the spin variables on the corners. It is shown directly that the row 
■ transfer matrices of the fused models satisfy special functional equations. 

Intertwiners between the fused A-D-E models are constructed by fusing 
the cells that intertwine the elementary face weights. As an example, 
we calculate explicitly the fused 2x2 face weights of the 3-state Potts 
model associated with the D4 diagram as well as the fused intertwiner 
cells for the A5-D4 intertwiner. Remarkably, this 2x2 fusion yields the 
face weights of both the Ising model and 3-state CSOS models. 



^ 1 Introduction 

The fusion procedure is very useful in studying two-dimensional solvable vertex and 
face models |], ||, [|. Essentially, fusion enables the construction of new solutions to 
the Yang-Baxter equations from a given fundamental solution. Among A-D-E lattice 
models [j|, || [7j , much effort has been focused on the fusion of the A models [||, H . By 
contrast, fusion of the D and E models has received no attention. The fusion procedure is 
important because it plays a key role in the solution of these lattice models. Specifically, 
it leads to solvable functional equations for the fusion hierarchy of commuting transfer 
matrices |)], Kj. Indeed, it has been argued [Tl[ that the fusion and inversion hierarchies 



of functional equations for the D and E models are exactly the same as those for the 
associated A model related to it by an intertwining relation []!□], I3J, [T3| . 



Here we extend the fusion procedure to all the critical A-D-E and the elliptic D 
lattice models. In particular, we establish the fusion and inversion hierarchies directly 
for the classical D, E and the elliptic D models. We also extend the construction of 
intertwiners to the fusion A-D-E models. In this paper, for simplicity, we focus on the 
classical A-D-E models although similar arguments apply for the affine and dilute A- 
D-E models. The paper is organized as follows. In the next section we define the critical 
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classical A-D-E lattice models and the elliptic D L models and modify the face weights 
by an appropriate gauge transformation. The modified face weights satisfy a group of 
special properties which ensure that they can be taken as the elementary blocks for fusion. 
In section 3 we give the procedure for constructing the fused A-D-E face weights. This 
is accomplished by introducing parities for the fusion projectors. In section 4 we derive 
directly the fusion hierarchies satisfied by the fused A-D-E row transfer matrices. The 
intertwiners between the fused A and the fused D or E models are presented in section 5. 
Also, in this section, we find the gauge transformation to obtain the symmetric fused face 
weights. In section 6, as an example, we give explicitly the fused D^ face weights and 
the fused cells that intertwine them with the fused A§ face weights. Finally, after a 
brief conclusion, we present in the appendices a comprehensive table of the adjacency 
diagrams for the classical A-D-E fusion models as well as the parities of the first four 
fusion levels of the E 6 model. 



2 Properties of the Face Weights 

The A-D-E lattice models |5], IB, |T7| are interaction-round-a-face or IRF models |18] that 
generalize the restricted solid-on-solid (RSOS) models of Andrews, Baxter and Forrester 
@]. At criticality, these models are given by solutions of the Yang-Baxter equation 



based on the Temperley-Lieb algebra and are associated with the classical and affine A— 
D-E Dynkin diagrams shown in Figure 1. States at adjacent sites of the square lattice 
must be adjacent on the Dynkin diagram. The face weights of faces not satisfying this 
adjacency condition for each pair of adjacent sites around a face vanish. 

In this paper we will restrict our attention to the classical A-D-E models. The face 
weights of the classical A-D-E models at criticality are given by f| 
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where u is the spectral parameter and A = n/h is the crossing parameter. Here 

L + l, for A L 

2L - 2, for D L (2.2) 
12,18,30, for E L = E 6 , 7fi 

is the Coxeter number and S a are the elements of the Perron-Frobenius eigenvector S of 
the adjacency matrix A with elements 
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otherwise. 



(2.3) 



In analogy to the classical A models, we modify the A-D-E face weights (|2.1| ) by a 
gauge transformation as follows 



g(d,c)g(c, b) 
g(d,a)g(a,b) 




where we set g(a, b) = g a g b with g a = S\/ A f^ 2 and 
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Figure 1: Dynkin diagrams of the classical and affine A-D-E Lie algebras 
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In this gauge, the modified face weights are given by 
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where we have introduced the symmetric sign symbol 



1, a = c 

1, (a, c) - 

1, (a,c) = 

— 1, otherwise. 



A _ I 1, (a, c) = (L - 1, L) or (L, L - 1) for £> L 
6a,c 6c,a / a ] 1, (a,c) = (L- 4, L- 2) or (L - 2, L- 4) for E L 



The face weights ( |2.1| ) or ( |2.8j ) satisfy the Yang-Baxter equations 





/< u—v 



(2.10) 



where the solid circles indicate sums over the central spins. 

Each node a of the A-D-E Dynkin diagrams has a coordination number or valence 
val(a) = 1, 2, 3. Specifically, the valence val(a) = 2 except for the endpoints with val(a) = 
1 and branch points with val(a) = 3. In the modified gauge ( |2.4|) the face weights acquire 
the following properties: 
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A a)b A b<c val(6) > 1. 



Moreover, at u — —A, the face weights also satisfy the properties: 
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These properties are useful for constructing the fused face weights. However, to study 
the fusion hierarchy we also need the additional properties: 
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We will introduce the corresponding antisymmetric sums in Section 3. 

However the fusion procedure constructed in Section 3 is described by studying the 
classical ADE models. In fact it works also for the elliptic Dl models with the nonzero 
face weights Wn which are related to the face weights Wa of the elliptic A 2 l-3 models 
by orbifold duality [|l^, [2f| 



W D 



L 
L 

W, 



u 



D 



W, 



D 



2 L - 1 

3 L-2 
L-2 L 
L-3 L-2 

L-2 L-l 
L-l L-2 
L-2 L 
L L-2 
L-2 L-l 
L-l L-2 



u 



L-2 L-l 
L-3 L-2 



u 



(2.27) 



u 



u 



A 



U 



2 L-l 
1 L 



u 



(2.28) 



5 



w D 



u 



w, 



D 



w D 
\w A 



w D 
w A 



II 



D 



w D 
w D 
w D 



L-2 L 
L-l L-2 

L-2 L-l 
L L-2 
L-2 L-l 
L-l L-2 
L-l L-2 
L-2 L-l 

L L-2 
L-2 L 
L-l L-2 
L-2 L-l 
L-l L-2 
L-2 L 

L L-2 
L-2 L-l 
L-l L-2 
L-2 L-l 
L-2 L - 3 
L - 3 L-2 
L-2 L - 3 
L-3 L-A 
d c 



a 



u ) - -W A 



L-2 L-l 
L-l L 



a 



u)-W A 



L-l L-2 
L L-l 



u 



ti 



W D 

w A 



II 



u 



it. 



uj +W A 

--w/' 



W A 



L-l L-2 
L L-l 
L-2 L-3 
L-3 L-2 
L-2 L-3 
L-3 L-A 



u 



a 



a b 



u 



W A 



d c 
a b 



ii 



if d^ L-2, L-l, L. 



Here the nonzero face weights W A are given by 
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where h(u) = 8i(u)8^(u), w a = aX and 9\, 64 are the usual theta functions of nome p. 

We have the same properties as ( 2.11 )-( 2.16|) if in the gauge transformation ( p.4|) we 
set g a = h l J 4 fl^ 2 where 
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and f a is given by ( |2.6| ). With these changes, the fusion of the elliptic D L models proceeds 
as for the critical A-D-E models. 



3 Elementary Fusion 

The Temperley-Lieb A-D-E models are related to the six-vertex model and hence to 
the spin algebra su(2). The higher-spin representations of this algebra are obtained by 
taking tensor products of the fundamental representation. The analog of this process for 
the A-D-E face models is fusion. Starting with a fundamental A, D or E solution of 
the Yang-Baxter equations it is possible to obtain a hierarchy of "higher-spin" solutions 
by fusing blocks of faces together. The fused A models have been discussed by a number 
of authors [0, 0, |9|, |], [L(| . In this section, we extend the fusion procedure to the classical 
D, E models and the elliptic D models. We focus on the critical ADE models and the 
arguments apply for the elliptic D models by replacing all sin u functions with the elliptic 
functions h(u). 



3.1 Admissibility 

The adjacency matrices A^ of the level n fused models are determined by the su(2) 
fusion rules jnj truncated at level h — 2 

A (n) A W = A (n+l) + A (n-1)^ n = ^ 2, 3, • • • , A - 2 

A W = I, A {1) = A, A (n) = 0, n>h-2, (3.1) 



A (h-2) 



I, for D 2 l, E 7 and E 8 
Y, for A L , D 2L -i and E 6 



where I is the identity matrix, h is the Coxeter number and Y is the corresponding height 
reflection operator defined by 

Y a ,b = $a,r(b) (3-2) 

where 



r(b) = h-b for A L (3.3) 
r(6) = l 6 ~ b if f ? <6 forE 6 (3.4) 




r{b) = \2L-l if b = 2L - 2 for £> 2i -i (3-5) 



Here A^ 1 ' = A is the adjacency matrix for the elementary classical A-D-E model. As 
examples, we draw the adjacency diagrams describing the allowed or admissible states of 
adjacent sites of the fused D 7 and El models in Appendix A. In contrast to fusing the 
Al models, the elements of A^ can in general be nonnegative integers greater than one. 
In this case we distinguish the edges of the adjacency diagram joining two given sites 
by bond variables a, ft — 1, 2, ... If there is just one edge then the corresponding bond 
variable is a = 1. 
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3.2 One by two fusion 



We implement the elementary fusion of a one by two block of face weights. The properties 
of this elementary fusion then suffice to establish the fusion of general m x n blocks of 
face weights. Notice that in the level 2 fused D and E models, the occurrence of bond 
variables on the edges of the fused face weights only arises when both adjacent sites are 
branch points with valence val(a) = 3. 

Lemma 1 (Elementary Fusion) If(a,b) and (d,c) are admissible edges at fusion level 
two we define the 1x2 fused weights by 
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(3.6) 



where the sum over a' is over all possible spins (i.e. a normal sum with a = 1) if a and 
b are not both of valence 3. If a and b are both of valence 3, the sum is accomplished in 
two different ways by summing over sym(a') . Explicitly, for Dl (resp. El) we sum over 
L — 3 and L — 1 (resp. L — 4 and L) if the bond variable a = 1 and over L — 3 and L 
(resp. L — 2 and L) if the bond variable a = 2. Then it follows that: 

(i) The RHS is independent of d except for its dependence on the bond variable 



c = d = L 
c = d = L 
otherwise. 



2 and d 

3 and d 



L: 



2; 



E L 



(3.7) 



(ii) For all a,b,c,d we have Wul 



d f3 
a a 







0. 



Proof: To establish (i) it is enough to consider the case c = d, otherwise d is uniquely 
determined by the adjacency conditions. Setting v = X and c = e in the Yang-Baxter 
equation ( 2.10|) we have 





(3.8) 



L + X 



Figure 2: Elementary fusion of two faces. The cross denotes a symmetric sum labelled 
by a = 1, 2 as defined in lemma 1. The other spins are fixed. If val(c) = val(rf) = 3 we 
assume that d ^ L — 3 for Dl and d ^ L for El- For clarity both the spin d and the 
bond variable (3 are indicated. 
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If a = b, then take the special sum over a' in ( |3.8| ). Owing to fl2.13|) , the special 
summation over a' with each fixed d vanishes in the LHS. Therefore for any (a, b) we 
always have 
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for D L 
for E Tl 



(3.9) 
(3.10) 
(3.11) 
(3.12) 



These equations imply part (i) of the lemma. Part (ii) follows by ( |2.11| ) if d ^ a and by 
( |T3D if d = a. 

Lemma 1 gives the 1x2 fused face weights incorporating the level two fusion adjacency 
conditions. A bond variable a has been added between each pair (a, b) of adjacent spins 
to form edges with states (a, a, b). The adjacency condition for bond variables is that 
a — 1,2 ii a = b = L — 2 (resp. L — 3) for D L (resp. E L ) and otherwise the bond variable 
takes the fixed value a — 1. Similarly, the spin variables are constrained by \a — b\ = 0, 2 
and 2 < \a + b\ < 2L - 4 (resp. 2L - 2) or (a, &) = (£- 1, L) (resp. (a, &) = (£- 4, L)) 
for (resp. £x). Observing properties (|2.14|) - (|2.17|) we find that this adjacency is 
completely determined by the operator P(l, — A) with elements 

d 



P(l, 



\\d,c,b 
~ A )d,a,b 




(3.13) 



So it can be considered as the projector of level 2 fusion. 

3.3 Operator P(n,u) 

Let us define graphically 
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r \' l i u )a,b 1 ,b2,-,b 




(3.14) 
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Then the operator P(n, —nX) is the projector of level n + 1 fusion. 

For n = 1 it is the face weight of an elementary block. For n = 2 it produces the 1 
by 2 fusion presented in the last section. This follows from the properties (|2.14| )-( |2.17| ) 
and ( p,13| ) we have 



P(2,u 



,d,c',c,b 
' d,a,ai,b 



P(l,-X) 
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a,ai, 



w 



12 



d (3 c 
alb 

' d (3 

12 I r, 

2 

d (3 c 

1 b 

d (3 

a 2 

(3 c 

1 b 



u} + 
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if a=b=L — 3 for E L 

if a =b=L-2 for D L 
otherwise 



(3.15) 



where a' is determined by the adjacency condition A a>a i = A a i h = 1. 

We now study the operator P(n, —nX) for level n + 1 fusion. With the help of Yang- 
Baxter equation ( |2.10| ) we can show that this operator satisfies 




(3.16) 




(3.17) 



These properties will be useful in later sections. 

Using the YBE ( |2.10| ) and the relations ( |3.15| ) it is easy to see that any two adjacent 
faces with spectral parameters u + jX and u + (j — 1)A in ( p,14j ) can be considered as an 
instance of 1 by 2 fusion. So the properties Q3.9p-( |3TT2j ) imply 



P(n, u] 



(a,bi, 



>i-l,fri;' 



-U-,-,b) 
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0, if val(6j 



val(6. 



i+ij 



(3.18) 
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Pf-n „\(a,ai,-,L-2,L-3,L-2,-,b) _ p( s(a,ai,-,L-2,L-l,L-2,...,b) 
u )(a,b u -,-A-iMM+u-,-,b) - r V 1 ' U )(aM,-,-,bi-i,bi,h+i,-,-,b) 

+p(n, 'o::;:;::::;.^;-^; for d l c>,m 

P(r, n ,\( a > a l>-> L - 3 > L > L -3>-,b) _ tj( \(a,ai,-,L-3,L-4,L-3,-,b) 
a )(a,bi,-,-,bi-uh,bi+i,-,-,b) ~ - r ^ i ' ^(a,6i,-,6i_i,6 i ,6 <+1 ,.»,6) 

+p(n, «)SS;:;^^^:::;? 6) for e l (3.21) 

P(n ,,\( a ' a l> ,-, '°'-l' a *-l _1 ' a »+l' _ .,\(»>»l.-!«i-l.»i-l+l."i+l>-,6) 

^(o,6i,-,6 i _i,6 i) 6i+i.-,6) ~~ ^ ' U J(a,& 1 ,---,6 i -l,&iA+i,-",&) 

for Oi_i = a i+ i and val(aj_i) = 2 (3.22) 

Let p(a, b, n) represent the set of all allowed paths of n steps from a to b on the 
Dynkin diagrams excluding paths, such as in ( |3.18| ), which only give zero elements for 
the projector. Similarly, let Put) rje the number of paths in the set p(a,b,n). For 
convenience let p(a, b, n)i represent the i-th path in p(a, b, n) and p(a, b, n)ij be the j-th 
element of p(a, b, n)j. So we can rewrite the elements of the projector P(n — 1, u) to be 

The operator P(n— 1, u) is a square matrix and can be written in block diagonal form. 
By the properties CTHCT we may have |P(n - 1,«)^ 4 J = |P(n - 

or P(n - 1, = P ( n - i.^jfcSl + P (™ - 1 > u )p(oAn)l for an y P ath PfoM)* 
and suitable j and j. If so we treat the paths p(a,b,n)i and p(a,b,n)j as dependent 
paths. Otherwise the paths p(a, b, n)i and p(a, b, n)j are independent. Suppose there are 
ra|"jj independent equations deriving from the properties ( |3.20| )-( p.22| ), then there are 

A|"L = -P( " L — m (ai) independent paths in p(a, b, n) where A^ b s is precisely the element 
of the fused adjacency matrices given in (|3.1|). We denote these independent paths by 
a(a, b, n), a — 1, 2, • • • , There are many ways to choose the independent paths but 

they all lead to equivalent fused models. The remaining paths should satisfy 

P(n ~ 1, = £ 0£t) ^ " 1. * = 1, 2, • ■ ■ , m[J 6) . (3.23) 

a=l 

p{n - 1, u)^±:;Z::i; b b = o, « > & - 2 (3.24) 

The value of (ftf^bn) * s zero ^ ^ ne P a ^ n P( a > n )i * s independent of the path a (a, 6, n) and 
is +1 or —1 otherwise. According to ( |3.23| ) we can divide p(a, b, n) into A/™^ independent 
sets defined by 

p{n, a, a, b) = {{p{a, 6, n)^^ ^ 0}, a = 1, 2, ■ ■ • , A$ by (3.25) 



The first path in p(n,a,a,b) is a(a, 6, n), the i-th path is denoted by p(n,a,a,b)i and 

(a,fe,n) 



p(n, a, a, denotes the j-ih element of the path p(n, a, a, We call 0^" n ) the parity 
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of the path p(a, b, n)i relative to the independent path a (a, b, n). By ( |3.16| ) it is obvious 
that 



/.(<*.<*) 

(a,b,n) 

(,(*.«) 

\a,b,n) 



\a,b,n) 
(b,a,n)- 



(3.26) 

(3.27) 



Equation ( |3.24j ) holds because all paths in p(a, b, n) with n > h — 2 are related by ( |3.22j ) to 
P(n ~ ^4^;:^^*^ = with vaLC^O = val(6 i+1 ) = 1. As an example, 
we give explicitly the parities of the first four fusion level of the E§ model in Appendix 
B. 



From ( |3.2U| )-( |3.22| ) it follows that the maximum number of terms on the right hand 



side of (£23D is two. Let us set fg = P(n - 1, and t) = P(n - 1, uf^. Then 

in general, by (|3.22|) , we can divide the submatrix P{n — 1,^)^'^™) of P(n — l,u) into 
columns 
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/ 3 
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3 (n) 
(a, b) 
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and t? k (1 < k < Puj>)) can t> e expressed as 



.a , , 0',/3) ,/3 
\a,b,n) l k ~+~ ( P( a ,b,n) l k 



by ( |3.20|) -( j3.21|) . For the models only the first group appears and = 1. For the 



£>l and Ex models the second group is related to first group. It is easy to see that 
detP(n-l,u)j|";J;^ = and det P(n - 1, u) = 

This means that the matrix P(n — 1, it) or P(n — 1, w)p(a tlj * s reducible. The irreducible 

operator p(n — 1, it, a, 6) is obtained from the reducible one P{n — 1, uf^ a \ '™j by picking 
the independent elements as follows 



p(n - 1,0,6) 



where t^j = P(n — 1, — wA)^'^j. So ( |3.23|) can be written as 



(«,{>) \ 

l 

An) 
A (a,b) 



\ A (a,b) 



«(«) 
>,*>) 



l(n) 

,(») 



(3.28) 



(a,b) 



P(n - 1, -nA)^ = £ tfg&j p(n - 1, a, 



(3.29) 
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Finally, using ( |3.29| ), the operator ( |3.14| ) can be factorized as 



d(") 

(61,6) 



p(a, a„,n)ff 



p(a, a n , n)j jn 



a=l 



n) 



u 




u+(n-2)A 




1 p{b±,b, n)i2 


p(6i 





(3.30) 



This result implies that the fusion can be carried out if the operator p(n — l,bi,a) is 
invertible. The existence of the inverse operator p(n — 1, b±, a) -1 is shown in Section 5.2. 



3.4 General Fusion 

Let m and n be positive integers and define 



I[w lxn 

k=l 
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(d 


f3 


c 
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U = " 


u 


V a 


a 


b 


/ a 


a 



r (rn) 



(d,a) 

j=l a2,---,a m 



p(a, d, m)j,k+i afc+i v(b, c, m) k +i 
p(a,d,m)j t k a k u(b,c,m) k 



Here a = p(a,d,m)j i i, b = u(b,c,m)i, c = u(b,c,m) m+ i, 
(3 = a m+ i and the summation over a k ranges over a k — 1, 
1 x n fusion in turn is defined by 



d 



u—(m—k)X \ . 
p(a, d, m)j ) 



(3.31) 



An) 

' (p(a,d,m)j^ k ,u(b,c,m) k ) 



m+ l, a — Oil, 

The 



W lxn \ 



-(n) 



d (3 c 
a a b 
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fc=i 



u+(ib-l)A 



(3.32) 



The fused face weights ( |3.31| ) associated with a bond state (a, a, 6) are obtained by 
summing over the dependent paths within the set p(n,a,a,b). Similar ideas have been 
applied to the fusion of the Af£' models in fl9|, |27j. The resulting fused face weights 
depend on both the spin variables a, b, c, d and the bond variables a, /3, p, v. For the A L 
models these bond variables take only the value 1 whereas they take values for the 
adjacent spins a, b for the Dl and El models. For the At, models the fused face weights 
do not change at all if we change the paths p(m, b, 1, c)i to p(m, b, 1, c)j and p(n, d, 1, c)i 
to p(n, d,l,c)i. But, for the Dl and El models, we have the following lemma: 

Lemma 2 // the path /3(d, c, n) is replaced with its dependent path p(n, d, (3, c)j then the 
fused weight 



W„ 



I d j c 
p v 
\a a b 
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E ■ 



(d,c,n) 





( d 
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■) 




\ a 


a 


b 





(3.33) 
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q(d, a, m)j,2 



g(d,a,m) j!m _i 



q(d,a,rri)j, ri 



P(d, c,n) 2 



f}(d,c,n) n -i f3(d,c,n) r , 



u 

1 


h ■< 


u+(n-2)A 
P 1 


u + (n-l)A 
1 










u — (m — 2) A 


► 


► 1 

a + (n — m) A 


1 

u+(n— m+l) A 


U — (?71 — 1) A 


V < 


i 1 

u +(n— m— 1) A 


i 

u + (n— m) A 



c, m)2 



c, m) r , 



v(b, c, m) r , 



p(a,b,n)i > 2 



p(a,6,n)i„_i p(a,b,n) iin b 



Figure 3: Diagrammatic representation of the face weights of the m x n fused ADE 
models. Sites indicated with a solid circle are summed over all possible spin states. 



Similarly, if the path u(b, c, m) is replaced by its dependent path p(m, b, is, c)j then 



A™) 

E 

v>=\ 



Proof: Let us first consider 1 x n fusion 
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\ a 
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(3.34) 
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Ci-l 



Zv ^(a,b,n) 












u — A 








p(a,b,n)j 


i+1 p(a,6 







(3.35) 



From ( |3.30|) it follows that the indices (c i+ i, q, of the weight ( |3.35| ) satisfy the 
properties (|3.9|) - (|3.12[) (or ( [3.19| )- (|3.22| )). This means that some of the fused weights in 
( |3.35| ) are dependent. In total there are A " « independent paths in the set p(d,c,n). 
Choosing an independent path q = /3(d, c, n) we have the lxn fused face weight 

d (3 c 



W lxn \ 



a a b 



a 



pC«) d 

A*,b) 

V( a ,b,n) 


/3(d, c, n)2 


m 


c,n) n 


u 




u+(n-2) A 


u+(n-l)A 


i 

a 


p(a, b, n)iz 


p(a, b, n)^ 



(3.36) 



where a = 1, 2, • • • , ij^ and j3 = 1, 2, • • • , ^.("dy These represent the independent fused 
face weights. The others can be obtained from the independent weights via the relation 



W lxn 
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o(") 
(a,b) 



Z^ < P(a,b,n 



(a,b,n) 

c 
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p(n 




p{n,d,P,c)ip, 


u 




u+(n-2) A 


u+(n-l)A 


p(a 


b, n)# 


p(a,b,n)i fl 



An) 

E 



d ft' c 



u 



(3.37) 



where ^^cn) are ^ ne P ar hies of the path p(n,d, ft,c)j relative to the dependent paths 

f3'(d, c, n) with ft' = 1, ■ ■ ■ , ft, - ■ ■ , A$ c \ . These properties are exactly the same as ( p.29| ). 
Furthermore, we have the following push through property from (|3.37|) 
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p(a,6,n) j;i _ip(a,b,n) ^ p(a,fe,n) Jiif i 



A ( E ,/) 
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f3=l 



p("A")j,HI | (»i l i»)j,if | Wi»)j 1 itl 





D-A 


D 




9(e,/,n)fe,i-l l(ej,n) k i q(e,f,n) k m 




w— A 







(3.38) 



This relation and ( |3.37| ) imply fl3.33|) . Moreover, (|3.34j) follows from ( p.33| ) because of 
d c 



the symmetry W 7 ^ , 



u ) = W 



b c 
a d 



u 



By repeated use of (|2.10| ), and with the help of the Lemma 2, we obtain the following 
theorem: 



Theorem 1 For a triple of positive integers m,n, I, the fused face weights ( \3. 3\ ) satisfy 
the Yang-Baxter equation 



/ f Vi 9 

E E^xn P V2 

im,m,m) 9 \a a b 



(e p d 
m 7 
9 Vi c 



d 7 c 
u | W mx i | r/ 3 ft 
9 m b 



e 7/2 g 

U | W mx i | V 7/3 

f P a 



e p d 

V-U ) W mxn \ V 7/3 

/ 7/i g 



9 7/1 c 

7J-71 | W mxn | 7/3 /3 

a a b 



(3.39) 



By Lemma 1 the weights W mxn ( p.31|) have zeros independent of the spins a, b, c, <i 
and bond variables a, ft, \i, v. To remove these zeros we replace the (M, N) fused weight 
by 
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n-2m-l • \ 

n n — 



(3.40) 
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(d 




c 




By construction it is obvious that W mxn 


n 




V 


■) 




y a 


a 


b 





u vanishes unless 



A^l ^ and a = 1, 2, 



.4 



(n) 



o,6 

19... 4^ 



4^0 and /3 
Ag/O and a* = 1,2 
Ag^O and z/=l,2,---,4J 

where the fused adjacency matrices are given by ( |3.1| ). In particular, 

if n = h — 1 or m — h — 1 





fd 
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•) 




\ a 
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(3.41) 



(3.42) 



4 Row Transfer Matrix Fusion Hierarchy 

Suppose that a(ot) and b((3) are allowed spin (bond) configurations of two consecutive 
rows of a lattice with N columns and periodic boundary conditions. The elements of the 
fused row transfer matrices T^ m,n ^(u) of the fused A-D-E models are given by 



N 



a i+i Vj+i bj+i 

(a,a\T m ' n (u)\b,(3) = ]JJ2 W mxn\ a, (3 3 

2 =1 {Vj} \ dj T]j bj 











Off! 




u 

) 


a 3 


u 









Pi 



(4.1) 



where a^+i = a>i, &v+i = &i and = rji. Specifically, the Yang-Baxter equations 

( p. 39 ) imply the commutation relations 



[T m ' n (u),T m ' n '(v)] = 0. 



(4.2) 



Thus if m is held fixed we obtain a hierarchy of commuting families of transfer matrices. 
These transfer matrices satisfy the following remarkable functional equations: 



Theorem 2 (Fusion Hierarchy) Let us define 



{u + kX), T 



m,0 



/-ll, fn=[s 



mlN 



and 



Then 



3=0 



sin A 



rpm,nrpm,l 
A A n 



fmrpm,n— 1 , rpm,n+l 
Jn 1 "r in-1 A 



where the hierarchy closes at fusion level h — 1 with 

r£P,h-l _ Q 



(4.3) 
(4.4) 
(4.5) 
(4.6) 
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Theorem 3 (TBA Hierarchy) If we further define 

mm,n+lmm,n-l 

^m,n 1 

— fm fm 

Then the thermodynamic Bethe ansatz equations 



(i+C' n+i )(i+tr n_1 ; 



hold where 



l m,0 



,m,h—2 
C 



0. 



(4.7) 

(4.8) 
(4.9) 



The main purpose of this section is to prove these theorems. Clearly, the functional 
equations for the Dl and El models are the same as those for the Al models. In the 
Al case the fusion hierarchy of functional equations was obtained by Bazhanov and 
Reshetikhin ||. Although intertwiners can be constructed [15[] between the row transfer 
matrices of the D or E models and an associated A model, these intertwiners do not 
relate all eigenvalues. Rather, only a subset of common eigenvalues are intertwined. As a 
consequence, the functional relations of the Dl and El models cannot be obtained from 
those of the Al models using intertwiners alone. Instead it is necessary to prove these 
functional equations directly for the Dl and El models as is done here. 

In Section 3 we described fusion of the A-D-E models corresponding to the symmetric 
representation of the tensor products of n elementary blocks. To prove the theorems we 
need the fusion procedure corresponding to antisymmetric representations. We therefore 
now describe the antisymmetric fusion of the tensor product of 2 elementary blocks. The 
symmetric and antisymmetric fusion procedures are orthogonal to each other in the sense 
that 



E E 

cGantisy e£sy a 



i+X 



(4.10) 



From ( p.l0| )-( |3.12| ) we can indeed see that (|4.10 ) holds where the antisymmetric sum is 
defined by 

L L-2 l 

d = L for Dl 
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cdantisy o, 
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otherwise. 



(4.11) 
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Furthermore ( |4. 1U ) implies that 



E 

c£antisy a 



i + A 



unless a = 6. 



(4.12) 



Hence, for the Dl models, we can construct the antisymmetric fusion by 

2 g 2 



& a 6 
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- E 

cdantisy 1 
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6,c 



-2 
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a = 1,6 = 2 



a = L,b = L-2 



otherwise. 



Similarly, the antisymmetric fusion for the El models is given by 

2 c 2 
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a = 1,6 = 2 
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By direct calculation we have 



a = L-3, 6= L-2 



a = L-4,6 = L-3 



a = L, 6 = L — 3 



otherwise. 
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1 < 
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u + A 



Oh 



i+A 



^a.c^l^— 1 



(4.13) 



(4.14) 



(4.15) 



Proof of Theorem 2: For simplicity we prove the functional equations only for the case 
of m = 1. The general case can be proved similarly. Representing Tq^T^' 1 graphically 

as 
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A {n) 
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p(a,b',n) iy 2 p(a, 6',n)i,3 p(a, f)', 




b',n)2 p 


(a, 6',n)3 
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« + A 











and inserting ( p.21[) we obtain the sum of two terms: 
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(2 cos A)" 1 
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(2 cos A)" 1 



p(a,b',n) it 2 p(a,b',n)i, 3 
fj.(a,b',n)2 fi(a,b' ,n):j 



p(a,b',n)i^s 
At(a, 6', n) n 



(4.16) 



(4.17) 



(4.18) 
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But now, by (|2.12|) , the second term ( [4.18| ) vanishes unless p{a,b', 



n 



i, n 



b. In this 



case we can choose an independent path with fi(a,b',n) n = b so ( }4.18|) becomes 
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Using (pl^)-pl^), (ggg)-(gjg) and (p^|)-(p^), this can be reduced to 
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(4.19) 



(4.20) 



By virtue of (|4.15|) this gives the first term /^Tq'™ 1 in the fusion hierarchy. From 



the push through property (|3.38| ) and ( |3.15| ) of the 1x2 fusion we can see that the path 
of 3 steps (fi(a,b f ,n) n ,b' ,b) in (|4.17|) satisfies the properties ( |3.9|) -( |3~7L2"|) . This together 
with the push through property (|3.38|) ensures that the path of n + 1 steps from a to b 
to satisfies ( |3.19|) -( [3.22| ). Applying the push through property ( [3.381 ) to the n + 1 blocks 



we obtain the level n + 1 fusion transfer matrix given by the second term /^jT, 



l,n+l 

o 



Proof of Theorem 3: Following Kliimper and Pearce |10| the functional equations 

rpm,7irpm,n jm y m J _|_ r-j-im,n+lrpm,n— 1 ^ 



are derived by substituting the fusion hierarchy ( |4.5| ) into the identity 



rpm.ra^rj-iri 



m,n— lrrim,l\ 



^ rji m , n rjrm , 1 ^ rji m , n — 1 



(4.22) 
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This then yields 
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(4.23) 



The functional equations ( 
dynamic Bethe ansatz p0, [21] 



q) are identical in form to the equations of the thermo- 



23] . The fusion hierarchy for the Al has been solved 



T0|| for the finite-size corrections and hence the central charges, scaling dimensions and 



critical exponents. A similar analysis can be carried out for the Dl and El models. 

The functional equations of the elliptic Dl models can be obtained by straightfor- 
wardly replacing the sinw functions with the elliptic functions h(u) in Theorems 2 and 
3. The functional equations of the elliptic A model are given in ||. Here we have shown 
that the functional equations of the elliptic Dl model are identical in form to those of 
the elliptic A 2 l+z model. 

5 Intertwiners and Symmetric Fused Weights 



Here we extend the A-D-E intertwiners constructed in 151 to the fused A-D-E models. 



We build symmetric fused face weights and generalize the intertwining relation to apply 
directly to the symmetric face weights. We also construct the intertwiners between the 
row transfer matrices of the fused A-D-E models. 



5.1 Intertwiners 

Let A and G be adjacency matrices of an A and a D or E model respectively. These are 
square matrices with nonnegative integer elements. Then the adjacency matrix C is said 
to intertwine A and G if 

AC = CG. (5.1) 

In general C is a rectangular matrix with nonnegative integer elements. Similarly, there 
is an intertwining relation between the symmetric face weights W A of A model and the 



symmetric face weights W of the D or E models if WE 




b b' 



where 




d p r 



b V 



(5.2) 



(5.3) 
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is a family of cells labelled by four bond variables. Here the cells vanish unless the 
spins d, a are adjacent sites of A, the spins c, b are adjacent sites of G and the spins 
a, b and d, c are adjacent sites of the intertwining graph C . The bond variables C\(c2) = 
1, 2, • ■ • , C a ,b{Cd,c)i v — 1, 2, • • • , G C) b- These cells satisfy two unitarity conditions which 
can be written in the form 





(5.5) 



Using the adjacency intertwining relation ( |5.1|) and the fusion rules (|3.1|) it follows 
that the same intertwining relations hold between the fused adjacency matrices, that is, 



A (n) C = C Q{n)_ 



(5.6) 



We therefore expect to find fused cells that intertwine between the fused face weights. 
Let us perform the following gauge transformations for the cells 



\ St fa A 



sin 

\ qG fG 
a \ >~>b Ja 



(5.7) 



(5.8) 



Here we do not need the bond variables because they take the value 1 for unfused face 
weights. The transformed cells can be fused in the same way as the A models. The level 



n fusion of the transformed cells (|5.7|) is given by 



JL 



Cl 



c 2 



Cn-1 



d 



d 



(5.9) 



where the solid circles indicate a summation over all possible paths p(d, a, n) of the A 
model. The fused cell satisfies the same properties with respect to the path p(c, b, n) of 
the G model as does the operator P presented in ( |3.19|) -( [3.22|) . We can therefore restrict 
our attention to the independent paths from c to b of the G model with q = /i(c, b, n) i+ i. 
Applying the intertwining relation ( |5.2| ) to the m x n blocks we therefore obtain the 
intertwining relation between the fused weights W^ n and W^ n given by ( |3.31| ) 
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(5.10) 



Here summation is implied over each of the inner bond and spin variables. Alternatively, 
we can fuse the transformed cells (|5.8|) giving 



P («) 

(d,a) 

/ , S°(<f,a,Ti) 



Cl C 2 



d 



(5.11) 



The path p(c, b, n) in the fused cell is satisfies the same properties as the operator P of 
the A model presented in (|3.19|) and (|3.22|) . That is, it is independent of ci, C2, • • • , c n _i if 
the fused cell is nonzero. We thus have another intertwining relation for W^ n and W^ n 
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(5.12) 



where again summations are implied over inner bond and spin variables. From the 
unitarity conditions (|5 .4j) — (pT3|) it is easy to check the unitarity conditions for the fused 
cells 




(5.13) 




SbSd 
S c S a 



8d,d'8uy 



(5.14) 



where {C^Cl) is (C n ,C%) or (C%,C n ). The bond variables v = v' = 1 for ((%,(%) = 
(C n , Cj) because there is only one independent path between two spins of the fused A 
models. In such cases we discard the bond variable between adjacent spins. 

The fused cells are given by both ( |5.9|) and (|5.11|) . They give the intertwining relations 
( |5.10 ) and ( |5.12| ) respectively, either of which can be taken as the intertwining relation 
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between the fused face weights of the A and D or E models. However, the fused cells 
( |5.9| ) and ( |5.11| ) are independent. Since we need both the fused cells and their conjugates, 
the fused weights of the D or E models cannot be obtained from those of the A model 



and the fused cells ( |5.9| ) alone 



5.2 Symmetric weights 



The fused face weights given by (|3.31| ) are not symmetric, that is 
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\c V 


b 





(5.15) 



To symmetrize the fused face weights we need to apply a gauge transformation. 



(n) 



Afa b \ matrix p(n — 1, u, a, b) is nonsingular. This can be shown using intertwiners. Specif- 



Although the operator P(n — l,u,a,b) does not have an inverse matrix, the AT' 



x 



l (a,b) 



ically, from (|5.2| ), (|5.9|) - (|5.11|) and the properties (|3.19| )- (|3.22|) we have an intertwining 
relation between the operators p A and the p G , 




(5.16) 



E 



b V 



Here we have expressed the operator pin — 1,— (1 — ^IV,^'*) graphically as a 
triangle with 



i=l 



f A (a,n,b) = J] V^i/Mo.Wi 



i=l 



From these equations, and with the help of ( 5.13 )-( 5.14|) , we can easily obtain 




a a 




V b 



b b V 



(5.17) 
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and the inverse of the operator p 



a' 


a' 


a a 
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/ " 2<5 ' 1 ' a) 
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f 2 0,^,6') 
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6' 6 6 



b V 



(5.18) 



As a result we have shown that pin — 1, — (n — 1)A, a, 6) is nonsingular 

6 6 6 6 



E 




AM IP Plp a \ 





u =E AJ\? f\ U = s * 




(5.19) 



We use the square root of this operator to build the symmetric face weights from the 
unsymmetric ones. 



Theorem 4 Define the i'™^ x A^\^ matrix 



G(a, b, n) = \Jf(ol, b, n)p(n — 1, —(n — 1)A, a, b)F(a, b, n), 
where F(a, b, n) is the diagonal matrix 

F{a, 6, n) = Diag [ f(a, 1, b), ■ ■ ■ , f(a, A$ b) , b) }. 
Then the symmetric weights 

G(d, a, n)^G(a, b, n) a ^ 
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(5.20) 



(5.21) 



(5.22) 
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(5.23) 



(5.24) 



Proof: The symmetry ( 5.24Q is implied by the symmetry of the elementary face 



4 c 6 

'i — 
— b 



(5.25) 
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The symmetry ( |5.23| ) follows from the equality 

d d £ 



(a,b) (d,a) 

E E 

a'=l fj,'=l 



.(m) An) 

(c,b) (d,c) 

E E 

u'=l (3>=l 




d d 






fJL 
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v 1 





f(d,(3,c)f(c,u,b) 
f(d,ti,a)f(a,a,b) 




which follows from the Yang-Baxter equation (|2.10 ) and 



(5.26) 
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(5.27) 



It should be noted that this gauge transformation is different from that used by Date 
et al || for the fused A models. Obviously, the symmetric fused weights W^ n and W^lxn 
satisfy the intertwining relations ( |5.10|) and ( 5.12 ) with the cells replaced by 



cub 
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G (d, a, n) 



G G (c,b,n) 
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G G (rf,a,n) M ^ 
V G A (a } b,n) 
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(5.29) 
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5.3 Row Transfer Matrix Intertwiners 



In this section we study intertwiners relating the row transfer matrices of the fused A- 
D-E models. 

It is easy to show [Fj| that the adjacency intertwining relation ( |5.6|) 

A (n) ~ G in) (5.30) 

is an equivalence relation among symmetric matrices. The existence of an intertwiner 
reflects a symmetry relating the two graphs associated with and G^ n \ In particular, 
we observe that the intertwining relation implies that 

[CC T , = [C T C, G {n) ] = (5.31) 

so that the symmetry operators CC T and C T C commute with A^ and respectively 
and their eigenvalues can be used as quantum numbers labelling the eigenvectors of A^ 
and GW. 

The above properties of intertwiners at the adjacency matrix level carry over to those 
at the row transfer matrix level. Let us introduce cell row transfer matrices with fused 
cells C n 



f>1 b 2 



a|C (n) |6), (a\Cf n) \b) 



N+l 



a N+l 



(5.32) 



where a and b are allowed row configurations of the graphs A^ and G^ 1 ' with periodic 
boundary conditions cln+i = a\ and b^+i — b\. In general, the row intertwiner is a 
rectangular matrix. Using the cell intertwiner relations (|5.10|) and ( 5.12|) it can be shown 
that 

A {n \u)C {n) = C {n) g^(u) or A {n \u) C £> G in) (u), (5.33) 



where 



i(n,m) 



U, 



and 



i(n,m) 
G 



u) are the row transfer matrices of two fused 



models. This intertwining relation can be pictured as follows 











On 





















N+l 



JV+1 - a l 

f _ f 

N+l ~ a l 



a N+l = a l 

where a solid circle indicates a summation over the corresponding spin. 



(5.34) 



27 



This intertwining relation is (i) reflexive, (ii) symmetric and (iii) transitive, that is, 
the intertwining relation is again an equivalence relation 

(i) A {n) {u) I A (n \u) 

(ii) AW(u) C ® gW(u) implies g^ n \u) C h A {n \u) (5.35) 
(iii) A {n \u) C £> B {n) (u) and B (n \u) C{ ^ G {n) (u) 



implies A {n \u) Q^(u). 



and moreover 



[C (n) Cf n) , = [Cj n) C {n) , <?<»>(«)] = 0. (5.36) 

Hence the symmetry operators C( n )C^ and C^C( n ) and the row transfer matrices «4.( n )(u) 
and C/( n )(u), respectively, have the same eigenvectors and can be simultaneously diago- 
nalised. The eigenvectors that are not annihilated by the symmetry operators give the 
eigenvalues that are intertwined and are common to A^ n \u) and Q^ n \u). Since 

C (n)Cf n) ~ Cf n) C(n) (5.37) 

it is precisely the nonzero eigenvalues of these symmetry operators that are in common. 

Let us now consider the A L -D( L+3 y 2 fused models, with L odd, and define the height 
reversal operators 71 a and 1Z D for these models by the elements 

N N 

{a\7Z A \b) = ]J <Wfe), (a\H D \b) = JJ S ajAbj ) (5.38) 

3=1 3=1 

where for the A models r(b) = h — b and for the D models 

(b, for6 = l,2,...,(L-l)/2 

r(b) = l(L + 3)/2, for b = (L + l)/2 (5.39) 
( (L + l)/2, for b = (L + 3)/2. 

These matrix operators implement the Z2 symmetry of the models. It is easy to show 
that the fused cell row transfer matrices satisfy 

c {n) cf n) = i + n A , cf n) e {n) = i + n D . (5.40) 

The row transfer matrices of the fused A and D models commute with the corresponding 
height reversal operators. An immediate consequence of this is that the eigenvalues of 
A(u)^ = T^' n \u) and V(u)^ = T^' n \u) are in common if and only if the corre- 
sponding eigenvectors are even under the Z 2 symmetry. In particular, since the largest 
eigenvalue has an even eigenvector, the largest eigenvalue is in common and hence the 
intertwined models have the same central charge 

e= *LU-$l±*>). (5.41) 
n + 2\ h(h-n)J y J 
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Similarly, following Kliimper and Pearce ]TU| , it can be shown that the conformal weights 



(A r , s ,A 

of the excited states are given by [^(| . 

[ht-{h-n)sf 



rr 



Anh{h 



n) 



[s - l)(n - s + 1) 
2n(n + 2) 



(5.42) 



(5.43) 



where s and r label the rows and columns of the Kac table and Sq is the unique integer 
determined by 1 < s < n + 1 and Sq — 1 = ±(t — s) mod 2n. However, in contrast to 
the case of the A models, nondiagonal terms with A r )S ^ A r s occur for the D models. 

6 An example: D4 

In this section we find the 2x2 fused face weights of -D4 model and construct explicitly 
the intertwining relation between the A§ and D4 models. The -D4 model is an interesting 
example because it corresponds to the three-state Potts model. 

The adjacency matrices for the fused A 5 and -D 4 models are given by the fusion rules 
( ^.1[ ). The adjacency graphs are thus as shown in Figure 4. 

The adjacency graphs decompose into two groups for level 2 fusion. The symmetric 
2x2 fused face weights of the A 5 model are [H] 
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cos(w + A) cosu 
sin A 

sin(2w) 
a/2 sin A ' 

sin(ii + 2A) sin(w + A) 
sin 2 A 

cos(w + 2A) cos(w + A) 



sin A 



(6.1) 



cos A 



sin A 

sin(u + A) sinii 
sin A 

sm(u — A) sin(w — 2A) 
sin 2 A 



for group 1 and 



2x2 



2 2 
2 2 



u 



cos(m — A) cosii 
2 sin 2 A : 
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A,: - 1 2 . 3 . 4 . 5 



level 2 fusion A^ 2 ) 



level 3 fusion 



level 4 fusion A^ 4 ) 



9^ 



9-9 



4 3 2 5 



2 4 



Q 

3 



level 2 fusion Z>( 2 ) 






2' 



3' 




level 3 fusion D< 3 ) U 




level 4 fusion £>< 4 ) Q Q Q Q 

1' 2' 3' 3' 



Figure 4: The adjacency graphs of the fused A 5 and D 4 models. 
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2 sin 2 A ' 

sin(-u — A) sinw 
2 sin 2 A ' 

cos(w — A) cos-u 
2 sin 2 A ' 

2 v / 3sin(w + 2A) 



sin 2 A 

for group 2. These weights satisfy the following symmetries 



W. 



2x2 



d c 
a b 



u)=W< 



2x2 



d a 
c b 



u\=W. 



2x2 



6 — d 6 — c 
6 — a 6 — b 



u 



(6.2) 



(6.3) 



g D (a, b) for the D paths such that 



To obtain the fused cells we find a gauge transformation g A (a, b) for the A paths and 



= 5> A (d,a) 



/M)" 1 . (6.4) 
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The unitary conditions then take the form given by ( |5.4j ) and ( p.5| ). Dividing the fused 
cells into two groups, we find 



for the first group and 



o » ^ f 3, » .3 

aL_L 1 L_i 

xf ;d 

\ 5i— » il' 5i— ► il' 

3L-. B 7 3L-. Q' 

/ > 3 3 | ^ 

sLX 3LJ3' 



\ 4 



Or 

2 1— > — b' 2 1 — • — b' 

d ; d 
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C — t' 4 t_! — \ 2 i 



_L -1 1 

V5 I 1 1 



-1 1 
1 1 



-M 1 1 

^2 1 -1 



V2 I 1 1 



(6.5) 



(6.6) 



for the second group. 

These cells satisfy the unitary conditions (|5.4j ) and (|5.5|) . Hence, from the intertwining 
relation the 2x2 fused D face weights must be given JT3] in terms of the A face weights 
by 




(6.7) 



independent of the spin d. Inserting the fused weights and fused cells given above 
we find the fused face weights of the -D4 model. Explicitly, for the first group the nonzero 
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weights read 
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where a, b,c,d — 1', 3', 3 are distinct. For the second group the face weights are 



W 2 D X2 
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\ sin A 
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sin(2u — A) 



sin A 



cosu cos(m — A) 
sin A 



sinw sin(A — u) 
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(6.8) 
(6.9) 
(6.10) 



(6.11) 



(6.12) 



(6.13) 



(6.14) 



where the bond variables // ^ v = 3', 3 . It can be directly verified that these fused 
weights satisfy the Yang-Baxter equation. In fact, the first group gives precisely the face 
weights of the critical 3-state CSOS model 0. The second group gives the weights of 
the 8- vertex model at the Ising decoupling point. 

The 2x2 fused face weights of the model have been obtained here via the in- 
tertwining relation. However, precisely the same results are obtained by following the 
fusion procedure presented in Sections 3 and 5. Although we have concentrated in this 
article on fusion of the classical A-D-E models, the affine A-D-E and dilute A-D-E 
models |2^] can also be fused using these methods. Similarly, the methods are easily 
extended to fuse the elliptic off-critical D models. 
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Appendix: Fused A—D—E Adjacency Graphs 
Adjacency graphs of the fused models 




fusion level 1 and 9 fusion level 2 and 8 




fusion level 5 fusion level 10 
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A. 2 Adjacency graphs of the fused E 6 models 



1 2 3 4 5 

fusion level 1 (unfused) 
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fusion level 2 
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fusion level 3 
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fusion level 4 
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fusion level 8 
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fusion level 10 
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A. 3 Adjacency graphs of the fused E 7 models 




fusion level 7 and 9 f us i on i eve ] § 



1 2 3 4 5 6 7 
fusion level 16 
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A. 4 Adjacency graphs of the fused E 8 models 



2 3 4 5 
fusion level 1 and 27 



fusion level 2 and 26 






fusion level 5 and 23 



fusion level 6 and 22 
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fusion level 9 and 19 



fusion level 10 and 18 




fusion level 11 and 17 




fusion level 13 and 15 



1 2 3 4 5 6 7 8 
fusion level 28 
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B Appendix: Parities (j> of the E 6 

Fusion level 2: 
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Fusion level 4: 
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The others are given by 0^ 4 ) = 0(6,J 4 ) 
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